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SUMMARY 

A general method has been developed* using the methods of kinetic 
theory, whereby the surface temperatures of bodies can be calculated 
for steady flight at any speed in a rarefied gas. The particular 
solution was made for a flat plate j however, the calculations can be 
easily extended to bodies of arbitrary shape. 

It was found that the aerodynamic heating problem in the absence 
of solar radiation, that is, for the case of nocturnal flight, becomes 
of negligible importance at altitudes of 125 miles and higher and up 
to steady flight speeds of 36,000 feet per second. The effect of 
solar radiation, for the case of daytime flight, becomes increasingly 
important as the flight altitude is increased. At an altitude of 150 
miles and higher, solar radiation is the predominating factor that 
determines skin temperature. Owing to the strong effect of solar 
radiation on skin temperatures at high altitudes, the desirability 
of nocturnal flight is indicated in order to minimize skin temper- 
atures . 

In order to maintain low skin temperatures, it was found, that 
the angle of inclination of the body with respect to the flight path 
should be kept as small as possible. This may be accomplished in 
practice by designing the body to be finely tapered and by flying 
the body at small angles of attack. 

It is pointed out that skin temperatures may be reduced by 
insuring thermal contact between portions of the skin inclined at 
positive and negative angles with respect to the flight path. As 
much surface as possible should be inclined at negative angles. 
Practically, this may be accomplished by boattailing the body. 



2 


NACA TN No. 1682 


In the event that an internal skin— cooling system is employed, 
it is shown that the rate of internal cooling must he of the same 
order of- magnitude or greater than the rate at which heat is lost 
naturally by emitted radiation. If the cooling rate is below the 
natural radiation rate, cooling has little effect upon skin temper- 
atures. 

It Is shown that, in the case of a missile designed to fly over 
a wide range of altitudes and speeds, it is desirable to make the 
emLssivity of the skin as high as possible. This conclusion, how- 
ever, Is based upon a skin surface for which the emissivity. Is inde- 
pendent of the wave length of the emitted and absorbed radiant energy. 

A possible method of reducing surface temperatures is indicated 
by the decrease in skin temperature which accompanies a decrease in 
thermal accommodation coefficient. This phenomenon may be used to 
advantage if it is possible to decrease the accommodation coefficient 
by altering the surface characteristics of the skin. 


INTRODUCTION 

Missile or aircraft flight at extremely high altitudes introduces 
many new problems to the designer. One of the more important problems 
that arises is due to the skin temperatures that are attained at very 
high velocities. The severity of the problem is commonly exemplified 
by consideration of the disintegration by combustion of most of the 
meteorites that enter the earth's atmosphere at high speed. 

In reference 1, Tsien has summarized the meager amount of 
research that has been done in the field of rarefied gas kinetics as 
applied to aeronautics. It Is shown that the general field of aero- 
nautics may be divided into three major regimes, and that the division 
lines, although poorly defined at present due to lack of experimental 
data, are marked by the ratio of the mean free molecular path length 
to a characteristic body dimension. 

The three main regimes of aerodynamics, or gas kinetics, may be 
described very briefly as: 

1. The region of conventional aerodynamics in which the mean 

free path length is negligible in comparison with the body 
size 

The region of free-molecule phenomena in which the mean free 
path is large with respect to the body size 
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3. The so-called "slip flow" regime* intermediate "between the 
conventional aerodynamics and free-molecule regimes 

In conventional aerodynamic theory* the first "basic assumption 
that is made is that the fluid under consideration may be treated as 
a homogeneous continuous medium. This assumption certainly should 
become invalid at very high altitudes owing to the extreme rarefaction 
of the atmosphere. A measure of the atmospheric density is furnished 
by the average length of path (mean free molecular path) that is 
traveled by a molecule between successive collisions. The mean free 
molecular path has been computed by Johnson and Slack in reference 2 
for altitudes up to 150 miles. The result of their computations is 
shown in figure 1. The existence of a large mean free path should 
not be construed as an indication that the number of molecules per 
unit volume is small, in the ordinary sense of the word. For example* 
at a reduced density at which the mean free path is 10 feet* there 
are still 10 13 molecules contained in a cubic inch. 

The work of this report will be confined to a study of the heat- 
transfer processes occurrin g in the region of free molecule flow. 

The lower limit of the ratio of mean free path to a characteristic 
body dimension at which the derived theory becomes invalid cannot be 
stated at this time for lack of experimental evidence. The numerical 
computations* however* will be confined to altitudes of 75 miles and 
higher. The estimated mean free path at 75 miles altitude is 1 foot. 

The specific purpose of the present report was to develop a 
method whereby the surface temperatures of bodies traveling at any 
speed in a rarefied gas could be determined. It should be pointed 
out that, to the authors* knowledge, no experimental work has been 
done in high-speed rarefied gas streams. Consequently* the numerical 
values obtained should be regarded with some reservation until exper- 
imental evidence is accumulated concerning the validity of the assump- 
tions that are made in the analysis. 


SYMBOLS 


A. constant of integration in Maxwell’s equation 


C 

*\ 

c y > 

c 2 


total velocity* feet per second 


conponents of total velocity* feet per second 



MCA TN No. 1682 


base of natural logarithms, 2 . 7 l 83 »" 

energy of incident molecules, foot-pounds per square foot, 
second 

energy of molecules at plate temperature, foot-pounds per 
square foot, second 

energy of re-emitted molecules, foot-pounds per square foot, 
second 

energy of rotation, foot-pounds per molecule 

translational energy of molecules striking front of plate, 
foot-pounds per square foot, second 

translational energy of molecules striking rear of plate, 
foot-pounds per square foot, second 

vibrational energy, foot-pounds per molecule 


Maxwell’s distribution function, dimensionless 


gravitational acceleration, assumed constant at 32.2, feet 
per second squared 

number of degrees of freedom of motion 

normal solar radiation at top of atmosphere, 93-^ foot-pounds 
per square foot, second 

P4 

Boltzmann constant, 5.66 X 10 foot-pounds per degree F per 
molecule 

mass of one molecule, slugs per molecule 

molecular weight, pounds per pound-mole 

number of molecules striking front of plate, molecules per 
square foot per second 

number of molecules striking rear of plate, molecules per 
square foot per second 
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N 

No 

Q 

R 


Hi 


R r 


S 

T 


Ti 

T P 

T r 

U 
*x 
*7 > 
u z 

. 

Y 


m 


V x 

Vy 

Y z 


\ 

1 


number of molecules per unit volume of gas 

Avogadro’s number, 2.73 X 10 26 molecules per pound-mole 

heat removed from skin by an internal cooling system, foot- 
pounds per square foot, second 

i 

Universal gas constant, 15kk foot-pounds per pound-mole, °F 
absolute 

incident radiant energy, foot-pounds per square foot, 
second 

re-emitted radiant energy, foot-pounds per square foot, 
second 

surface area, square feet 
absolute temperature, °I' absolute 
temperature of incident molecules, °F absolute 
temperature of body skin, °F absolute 
temperature of re— emitted molecules, °F absolute 
mass velocity, feet per second 

components of mass velocity, feet per second 

thermal velocity, feet per second 

most probable molecular speed, feet per second 

components of thermal velocity, feet per second 
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X 

7 > 
z 

4 

a 

P 

€ 


space coordinates 

accommodation coefficients dimensionless 

reciprocal of most probable molecular speed, seconds per foot 
emissivity, dimensionless 


6 angle of incidence, degrees 

a Stefan-Boltzmann constant, 3.7^ X 10 —1 ° foot-pounds per 

square foot, second, °F 4 absolute 

2 P S —x 2 

$ (a) error function, -y— » / e dx, dimensionless 

a/ rt Jo 

X dimensionless group defined by equation (7) 

X* dimensionless group defined by equation (9) 

^ dimensionless group defined by equation (l4) 

dimensionless group defined by equation (19) 


ANALYSIS 

In order to determine the skin temperature of a body moving at 
constant speed through a rarefied gas, it is necessary, fundamentally, 
to mak e an energy balance on a section of the skin. Energy may be 
added to or subtracted from the skin by three distinct processes: 

(a) Molecular energy transport to and from the skin 

- (b) Eadiant energy transport to and from the skin 

(c) Energy- added to or removed from the skin by processes 
occurring within the body 

Thus, we may write a simple energy balance on a section of the skin 
as 


Ei + Hi = Er + Rr + Q 


( 1 ) 
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Each, term in this equation may he considered separately. 


Molecular Energy Transport to the Body 

It Is a matter of common knowledge that any body surrounded by 
a gas is constantly bombarded by molecules of the gas. It Is this 
molecular bombardment, in fact, that causes a pressure to be exerted 
on a body or on the walls of a container in which the gas is confined. 
It is evident that the molecules possess, by virtue of their mass and 
velocity, a certain kinetic energy. Thus, a body Immersed in a gas 
is constantly receiving and emitting energy from the bombardment of 
impinging molecules and the re-emission of the same molecules from 
the body surface. If it were considered that the molecules carried 
no energy except that due to translational velocity, then, we should, 
by examination of their incident and re— emitted velocities, be able 
to determine the net increase or decrease of the kinetic energy of 
the re— emitted gas and, hence, determine the molecular energy added 
to or subtracted from the body. In the ordinary case, however, the 
impacts of the molecules with the body become obscured by the col- 
lisions of the molecules among themselves and the problem becomes 
exceedingly complex if treated on a microscopic basis. On the other 
hand, if we consider a body immersed in a gas that is of low enough 
density so that collisions among molecules in the vicinity of the 
body are rare in comparison with collisions with the body, the 
problem becomes susceptible to analysis by the conventional methods 
of kinetic theory. 

If it is assumed that the gas through which the body is moving 
is so rarefied that the motion of the impinging molecules is essen- 
tially unaltered by collisions with reflected molecules, then it 
follows that the velocity distribution of the impinging molecules 
will be that of a gas in thermal equilibrium, the Maxwell distri- 
bution. The Maxwell distribution law for a gas at rest states that 
of a total of N molecules per unit volume, the fraction that has 
velocity components lying between the values, 

Yx to Vx + dYx, Y y- to Yy +dYy, Y z to Y z + <3Y Z is given by 


f(YxYyY z ) dYx®TydY z = Ae“P 2 (^x 2 +7y 2 +Y z 2 ) (2) 


where Y z , Yy, Y z are the random thermal velocities with respect to 
a coordinate system x,y,z fixed in space. 

To determine the energy transported to the body by the impinging 
molecules. It Is necessary first to determine the number of molecules 
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that strike the body. This is done as follows: 

Assume that the body is stationary and that the gaB has a 
macroscopic velocity U with components U x , U y, U z with respect 
to the fixed coordinates. Assume the body to be a flat plate and 
fix the plate at the origin of x, y, z, so that the x— axis is 
normal to the plate. The mass velocity U makes an single & with 
the plate. Figure 2 is a diagram of the coordinate system. 

The total, velocities may be written as the sum of- the thermal 
and mass velocities, thus, C x =U x +'Vx.> Cy=Uy+Vy, and C Z =TJ Z +V Z . 
Consider an elemental area dS of the plate. All of the molecules 
that strike dS during a time interval dt that have velocities 
lying in the range C x to C x +dC x , Cy to Cy4dCy, and C z to C z +dC z 
will lie in a cylinder having dS as a base and with a length C dt 
in the direction of C, where C= */C x 2 +C y 2 +C z 2 . The volume of the 
cylinder is C x dS dt. The number of molecules of this kind contained 
in the cylinder is then f'(C x ,Cy C Z )NC X dS dt dC x dCy dC z where 
f ' ( C x , Cy , C z ) is Maxwell's distribution function for the total veloc- 
ities. Then> to find the total number of molecules striking dS per 
second, we must integrate f '(C x ,Cy,C z )NC x dS dC x d Cy dC z between 
the limits C x =0 ton, Cy = -«to », and C z = — » to °° . Mole- 
cules having a negative velocity in the x-direction can n ot strike 
the front side of the plate. Then, 


n 



f'(C x ,Cy,C Z )C X 


dC x 


dCy dC z 


(3) 


The distribution function for the thermal velocities with a super- 
imposed mass velocity may be written as (cf., reference 3> P* 

f'(C x ,Cy,C 2 ) o Ae“^ 2 [( c srUz) 2+ ( c y -U y) 2+ ( c z -u z) 2 ] 


Thus for a unit area, 

n = NA P P P e~£ 2 [( c x^x) 2+ (Cy--Cry) 2 +(C z -U z ) 2 ] Cx dCx dCy dCz (4) 

d- do u— 00 d o 

This equation can be integrated; making use of the relation A * 
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the result is 


IT 


{ 


2 + 

2^ 9 e 


[ l+*(PUx> ] } 


(5) 


The factor 9 is related to the most -proba ble mo lecular speed 
as p =s while Y m is given, by V m = Also U* may be 

written as U sin 9. Equation (5) may then be rewritten as 



sin 2 9 

+ 


*/ir — sin 9 



This is the result obtained by Sanger in reference 4 B.nfl re— derived 
by Tsien in reference 1. A plot of the dimensionless group 
defined as 


( 6 ) 



To determine the number of molecules striking the rear side of 
the plate n*, equation (4) is again integrated except that the l imi ts 
of C x are changed to Cx = — 08 to 0. The result is 


n* = 


M m 

2*/* 


f-UJ 

l e 


sin 2 9 


/- U 

jy/n =— sin 9 
Vm 


1-4 isinS 


)] 


(8) 
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The group 


Y 


-\k) sln * 9 


X* =» < e 


is plotted in figure ■4-. 


/a sin 9 
v m 


1— $ ( =— sin 0 


(9) 


To determine the kinetic energy incident upon the plate due to 
translation, Sanger in reference 4, and Tsien in reference 1, multiply 
the mass of gas striking the plate nm by the sum of the kinetic 
energy per unit mass due to the mass motion of the gas TP f 2. and the 
average kinetic energy due to thermal translation which may be shown, 

3 HT 

for a static gas, to have the value -g-g -jj-. This procedure, however, 
is an oversimplification of the problem as introduction of the term 

3 RT 

assumes that the kinetic energy due to thermal velocity may- 
be added to the kinetic energy due to mass velocity. That this iB 
not permissible is explained by the fact that the total kinetic 
energy depends upon the square of the total velocity, hence the com- 
ponents of the kinetic energy due to the velocity components of the 
total velocity may not be added linearly. Symbolically, this may be 
expressed in simple form as . .. 

mV 2 mU 2 / m( U+Y ) 2 
2 2 f 2 

A more exact calculation of the incident translational energy can be 
made as follows: 


The kinetic energy per molecule contained in a cylinder placed 
on dS is mC 2 /2. The kinetic energy of those molecules in the 
cylinder with velocity components between C x to C x + d.C x , 

C j to Gy + dCy, C z to C z + dC z is 

j| NC 2 C x f'(Cx.,C 7 ,C z ) dS dt dC x dCy dC z (10), 

The total kinetic energy of the molecules striking the front side of 
the plate is obtained by integrating equation (10) over the limits 
c x = 0 to oo, Cy = — oo to oo, and C z = — 00 to 00 • Then, for a unit 

area, per unit time. 
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Et “ 2 N / / / c2c x f’(Cz,C y ,C z ) dC* dCy dC z (11) 

loi) 'ibinS - ^3 

fr eTer ' foK-r£&fdi? 4 t.- C-x^lfeMSSOJS^ (P-lse? Jft WIsSJ&w 

■tor.w .ae& ojclsj-a 0 ni . ovB-jn, spslq y’XBxriserer o/> saiaaoio Tgia?® 

asasBoiq-qa coiiai/pa add- ,9310! sJf <3 nia ! ’ 3 v\U ) 

Et = § HA f P P C 2 ( C z 2 +C v £ +C z 2 ) e - ^ 2 Ji C ^ -l:J; s) 2 + (Cy-D'y) 2 +(C 2 ;--D r z ) 2 3 

., *7- J-«o «X» Jo * ‘. T- -r •*• — ** » -t ■- vS 

! ■ J - ) \ il Q / 

dCx dCy dC z fl2) 

": : ‘r.?rr‘s Cflcoid-Blexteid- juahi oni erfj rto’t t ai>laiv. aoiStBTwqjso'.-. vjjliSS r«. a. 

The result of a tedious ^ ,is ’ S£9 ‘ '**** ««* .* 

. 01 ) v cni H' + — } f u ^ 't& 


Efc = nm 


I u! I 


Sf. + JL. 

, , 1 J V* PHj [l+®Ou*)J, e^x 2 ) 

,2 / 2p 2 ' 

1 + v^PUx [l+«(pux)] eP^x 2 

W V / f 

v '?« y t 

'■■' :-Tr^_± ) $-1 ! 6 ntB ~ A. X _ f 

given hy equation (6). Using the relations/ 


■j 2 kT' ^ ~ VZ’ = U sin 6 

r- i - / . 

■ ' . -T? ■ - ! 

i w ! 

and defining a dimensionless group 1 ty as 


U 

rnV 


( 13 ) 


» * ■* 
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When the mass velocity TJ is taken as zero, equation (15) reduces to 

Et = 2nkT ( 16 ) 

which is the usual result obtained for the translatory kinetic 
energy crossing an imaginary plane drawn in a static gas. When 
(u/V m ) sin 0 is large, the equation approaches 



A similar computation yields, for the incident translational energy 
on a unit area of the rear side of the platq,. 



The value of Et* , .of course, reduces to 

Et* = 2n*kT * 2nkT (20) 

for U equal to zero and approaches Et* = 0 for large values of 
(u/Vm) sin 0. The dimensionless groups ijc and are shown in 

figures 5 and 6, plotted as functions of (U/V m ) sin 0. 

If the gas is considered to be composed of hard— sphere molecules, 
then, for a monatomic gas, the total kinetic energy incident upon 
the plate is given by equations (15) and (l8) for the front and rear 
sides, respectively. However, if the gas is composed of diatomic 
dumbbell molecules, each molecule may carry an additional amount of 
energy due to rotation and vibration. By the principle of the equi— 
partition of energy, the maximum amount of energy carried by a 
molecule is kT where J is the number of degrees of freedom of 
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motion. A rigid dumbbell molecule may rotate about two mutually 
perpendicular axes, hence the maximum average kinetic energy of 
rotation per molecule is 


% = kT (21) 

It is conceivable that collisions between dumbbell molecules could 
cause the atoms composing the molecule to vibrate or oscillate along 
the line joining their centers. In this case, the energy of vibra- 
tion is composed of kinetic and potential energies and the vibra- 
tional energy is eq uall y distributed between the kinetic and potential 
forms. Thus, the total energy due to vibration is 

= 2(|kT) = kT (22) 

where the factor 2 accounts for the distribution of energy between 
the kinetic and potential forms. 

Some question may arise as to the correct procedure to follow 
when considering a gas composed of two or more constituents. It can, 
be shown, that the results that have been obtained for a homogeneous 
gas, may be extended to mixtures of gases by calculating n, Et, 

E r , and E v separately for each constituent. This is a result of 
the fact that, in a mixture of gases in equilibrium, each gas has the 
same velocity distribution that it would have if the other constit- 
uents were not present. 

The result obtained in equations (15), (16), (l8), and (20) 
shows that, in effect, both front and rear sides of the plate 
selectively collide with a group of molecules having an average 
thermal velocity in excess of the average thermal velocity of the 
molecules of the gas considered as a whole. This is true whether 
the plate is moving or stationary. As the plate velocity increases, 
the average thermal velocity of those molecules striking the front 
side approaches the average thermal velocity of the molecules 
of the gas considered as a whole, while the average thermal velocity 
of those molecules striking the rear side departs farther from the 
average thermal velocity of the molecules of the main mass of gas. 
Because of this selective interception, the plate is struck on both 
front and rear sides by molecules that possess higher— than— average 
translational energy associated with their higher— than-average 
thermal velocity. It is pertinent to inquire if this does not mean, 
also, that this nonaverage group of molecules carries nonaverage 
rotational and vibrational energies. It is pointed out in reference 
5 that Maxwell’s distribution law allows the assumption to be made 
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that rotational and translational energies are completely independent, 
so that a molecule with a high translational energy is just as likely 
to have a low as a high rotational energy. In the case of vibrational 
energies, it is shown that those molecules that have higher transla- 
tional velocities also carry higher vibrational energies. However, 
for reasons which will presently be discussed, the change in the 
vibrational energy is neglected in the energy balance equation, hence 
we need not be concerned with its magnitude in this simplified 
analysis. 


Molecular Energy Transport from the Body 

In order to determine the energy carried from the skin by 
molecules returning to the main mass of gas, it is necessary to 
consider the mechanism of energy transfer and molecular reflection 
at the surface of the body. A complete discussion is beyond the 
scope of this paper; the reader is referred to references 5 and 6 
for a detailed description of the phenomena involved. Nevertheless, 
the more important concepts are outlined below because an under- 
standing of them is necessary in succeeding computations. 

Two main types of reflection may occur when a molecule strikes 
a surface. The first and simpler type is specular reflection. In 
this case, a molecule leaves the surface with its normal velocity 
component reversed, retaining unchanged its tangential component. 

Any other type of rebound from a surface is known as diffuse reflec- 
tion. It is possible to devise many subtypes of diffuse reflection; 
however, only a limited number are physically probable or even 
thermodynamically permissible. Insofar as direction of reflection 
or re-emission is concerned, it can be shown that the probability of 
a molecule leaving the surface at a specified angle is proportional 
to the cosine of the angle with respect to the normal. All directions 
about the normal are equally probable, regardless of the direction 
of impingement. Several hypotheses may be advanced concerning the 
speed at which diffusely reflected molecules are emitted. The most 
physically likely assumption is that the speeds are grouped in a 
Maxwellian distribution corresponding to a temperature which Is not 
necessarily that of the surface. 

It appears logical that diffuse reflection would predominate as, 
in reality, an impinging molecule collides, not with a smooth hard 
surface, but- with an individual molecule or atom of the body material. 
It would seem likely that collisions of this type would obliterate, 
insofar as directions are concerned, any prior directional history 
of motion, especially if more than one collision occurred at the 
surface . 
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As yet, the temperature at vhich the molecules leave the 
surface has not been specified. In general, the molecules that 
leave the body will not possess an average velocity corresponding 
to the plate temperature . In the case with which we are concerned, 
namely, that of a cooled plate, the molecules will leave at some 
temperature which is higher than the plate surface temperature . 
This inequality of temperature (or energy levels) is defined by 
introduction of the accommodation coefficient, defined as 


a = 


e i -e p 


(23) 


The accommodation coefficient may be considered as a term that accounts 
for the inability of the impinging molecules to adjust themselves to 
the plate temperature dur in g the time they are in contact with the 
plate. As pointed out by Loeb, in reference 5, it appears that the 
action of surface roughness is different for directional effects and 
energy or temperature effects. Apparently, one or two collisions 
with the molecules of the plate is sufficient to erase the directional 
history of a molecule; however, many collisions are required to attain 
temperature equilibrium with the plate. In the definition of the 
accommodation coefficient given by equation (23), the effects of 
specular and diffuse reflection on the energy exchange are included 
in the value of a. 


The value of the accommodation coefficent has been measured 
for several static gases in contact with many different materials. 
Values reported range from less than 0.1 to 1.0, as shown in 
references J to 12. Reference 13 presents a method for correlating 
the effect of pressure upon accommodation coefficients. Apparently, 
no measurements have been made of the accommodation coefficient in 
a high-speed gas stream. It would appear probable, however, that the 
accommodation coefficient would depend upon both the speed of the 
gas stream and the angle of the surface with respect to the stream. 

The accommodation coefficient, if defined rigorously, would 
have different values for translational, rotational, and vibrational 
energies. However, a series of experiments performed by Knudsen 
and described in reference 5 showed that the accommodation coeffi- 
cients were numerically identical for rotational and translational 
energies. The vibrational components of the internal energy require 
much more time to adjust to new values than do the rotational and 
translational components. A discussion of the time of adjustment of 
vibrational energy for gases may be found in reference 14. In 
succeeding calculations, it will be assumed that the vibrational 
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component of- the internal energy is unchanged during the short time 
interval that a molecule is in contact with the body. By the same 
line of reasoning, dissociation of-polyatomic molecules is neglected, 
since dissociation is a function of the intensity of the molecular 
vibrations. These assunptions are serious in that - vibrational energy 
and the energy absorbed during dissociation would, if present', have 
a markecL effect on the plate temperature. If the length of time a 
molecule is adsorbed on the plate is equal to or greater than the 
time required for fully developed vibration to be ezcited, then 
dissociation would probably occur at the high plate temperatures 
calculated for high-speed flight. A more complete analysis awaits 
experimental data on adsorption times and accommodation coefficients 
in a flowing gas. 

To calculate the energy abstracted from the plate due to the 
return of molecules to the stream of gas, use is made- of equation 
(23) rewritten as 


E^ = { 1— ct) 4* ctEp (2^-) 

If Ep is the energy of the molecules if- they were at plate temper- 
ature, then equation (25) can be written for a monatomic gas that 
obeys the perfect - gas laws as 

Er = Et (l-tt) + \ ahkT p (25) 

and for a diatomic gas that obeys the perfect— gas laws, ■ neglecting 
vibrational energy changes, as 

E r = (E t + nkTi) (l-a) + \ ankT p (26) 

Radiant— Energy Exchange With the Body 

Any body is constantly receiving and emitting radiant energy. 

In the steady state, a balance is soon reached where the emitted 
energy is equal to the incident energy; the body temperature adjusts 
itself to a value such that this equilibrium is maintained. 

In the present case, the body emits radiation at a rate given 

by 


R r = eaTp 4 (27) 

The body may absorb radiation by several processes. Probably the 
most important— is that due to direct solar radiation. The mean 
value of the solar radiation at the top of the earth’s atmosphere, 
the solar constant J Q is shown in reference 15 to have a maximum 

1 
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value of 93 foot-pounds per second, per square foot for a plate 
normal to the incident radiation. The body may also absorb radia- 
tion from constituents of the atmosphere and from the earth’s surface. 
Most of the atmospheric gases do not emit appreciable radiation; 
however, water vapor, ozone and carbon dioxide emit strongly at 
certain wavelengths. In the interest of simplicity, radiation from 
the atmospheric constituents and from the earth’s surface will be 
neglected and, where incident radiation is considered, it will be 
assumed to be due solely to solar radiation, having its previously 
mentioned maxi mum value of 93 *^ foot-pounds per second, per square 
foot. The amount of this energy that is absorbed is then 

Hi = eJ 0 (28) 

provided that gray body absorption is assumed; that is, the emissivity 
is constant for all wavelengths and temperatures. 


Energy Exchanges Occurring Within the Body 

It is conceivable that, in a missile designed for sustained 
flight at very high speeds, a skin-cooling system could be used to 
maintain structure temperatures within specified limits. Disregard- 
ing a discussion of the possible cooling systems, we may simply 
assume that the cooling system removes heat at a rate Q foot-pounds 
per second, per square foot. 


CALCDIATIOES 

Having outlined the processes by which energy is exchanged with 
the body, it is now possible to calculate the temperature of a plate 
moving at any speed through a rarefied atmosphere. 


Method of Calculation 

We may substitute equations (25) or (26) , (27) , and ( 28 ) in 
equation (l) and obtain, for a monatomic gas, where E^ = Et 

oEt - ankTp + eJ Q - ecrT^ 4 - Q = 0 (29) 


and for a diatomic gas, where Ei = Et + h% 
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a (Et + nkTj) — j| ahkTp + eJ Q — eoTp 4 — Q = 0 (30) 


In equations (29) and (30), the value of E-fc may be t alien from 
equation (15) and the value of - n from equation (6) for the. front 
side of the plate. For the rear side of the plate, equations (l8) 
and (8) are substituted for values of E-^' and n', respectively. 

For a mixture of monatomic anri diatomic gases, the incident 
energies of the components of the gas are added together, and the 
re— emitted energies of the components of the gas are added together, 
so that the energy equation assumes the form 

EE^ — EEp + eJ Q — effTp 4 —0 = 0 (31) 

To determine the temperature of a plate with no cooling except- 
that due to radiation, the factor 0 in equations ( 30 ) and (31) is 
set equal to zero and the equations are solved for values of Tp. 

To determine the effect of internal cooling on the plate 
temperature, arbitrary values of Q are substituted into equations 
(30) and (31) and the equations are again solved for values of Tp. 


Scope of Calculations 

Calculations were made ter determine the temperature of a plate 
traveling at velocities of from zero to 36,000 feet per second at 
altitudes of from 75 to 190 miles. A speed of 36,000 feet -per 
second is approximately the minimum speed required to escape from 
the earth at these altitudes. The plate, in the first case, was 
assumed to be insulated between the-front and rear surfaces so that 
the two surfaces could be treated independently. In the second case, 
the front and rear sides were assumed to be in perfect thermal 
contact; that is, the temperature of- the front. side of the plate was 
equal to that of the rear side of the plate.- In these two Bets of 
calculations, it was assumed that solar radiation was absent. 

The above calculations were then repeated, for 6ne altitude, 
and single of incidence (or body angle) with arbitrarily assumed 
values of internal cooling and with solar radiation absent. 

A fourth set of calculations was made for several altitudes, 
considering that the plate was uncooled and that solar radiation, 
was present. 
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In the first four set? of calculations the values of the 
accommodation coefficient said the emissivity were assumed to he 
constant at unity. This assumption simplified the computations 
and, at the same time, reasonably approximated physical conditions 
so that the trends of skin temperature shovn should be correct in 
order of magnitude. 

A fifth and sixth set of calculations were made to show the 
effect of emissivity and accommodation coefficient on skin temper- 
ature at several altitudes and flight speeds. 

The atmospheric properties used in the calculation were taken 
from references 2 and 1 6. Reference 16 was used to obtain the 
properties of the atmosphere at 75 miles; whereas reference 2 was 
used for the higher altitudes. There is some difference in the 
assumed atmospheric compositions between the tables of references 2 
and 16. Warfield, in reference 16, assumed that the composition by 
weight at all altitudes was constant at its sea— level value, except 
for dissociation of oxygen. The amount of dissociation was assumed 
to increase with altitude, complete dissociation occurring above 
altitudes of 75 miles. On the other hand, Jo hns on and Slack assumed 
sea— level composition at n.1 1 altitudes, neglecting dissociation 
entirely. In order to compare values of skin temperature at the 
several altitudes for which computations were made, it was arbi- 
trarily decided to use the values of temperature and density given 
by Johnson and Slack and to assume complete dissociation of the 
assumed constant sea— level percentage of oxygen. The following 
table lists the values of atmospheric density and temperature that 
were used in the computations. 


Altitude 

(miles) 

Temperature 
(°F absolute) 

Density 
( lb/ ft 3 ) 

75 

680 

2.55xl0 -8 

100 

1025 

1.05x10“ 9 

125 

1385 

1.00x10- 10 

150 

1735 

1.75x10-11 

190 

1990 

1. 61x10“ 1 4 
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Discussion 

Incident Translational Energy- 

In the computation of the incident kinetic energy due to 
molecular translation it was shown that both the number of molecules 
striking unit area of the body and the resultant kinetic energy inci- 
dent upon the body were functions of the term (u/Vm) sin 0. An 
examination of the graphs of the dimensionless groups X and X’ , 
ijr and tS which govern the number of molecules striking a section 
of the body and the magnitude of the incident translational energy, 
reveal several interesting features. With respect to the number of 
molecules that strike a body on the front side, it can be seen from 
figure 3 that this number increases practically linearly with flight 
speed for a given plate or body angle. On the other hand, as shown 
in figure 4, the group X' that determines the number of molecules 
that strike the rear side of a body, decreases rapidly with flight 
speed and is practically zero for values of (u/V m ) sin 8 greater 
than 2. 

From equations (15) and (l8), it may be seen that the total 
incident translatory kinetic energy term is composed of two parts: 
the first due to mass motion of the gas relative to the body and -the 
second due to the apparent translatory molecular motions. In the 
case of the front side of the body, the apparent molecular energy 
term can be seen, from figure 5.> to increase rapidly at first with 
increasing values of (u/Vm) sin 8 and then to approach a constant 
value of 5/2 at values of (u/V m ) sin & greater than 2. On the 
rear side of the same body, the molecular energy term continues to 
decrease rapidly as (U/v m ) sin 0 increases. 

A simple physical interpretation of the trends noted above can 
be given. First, it is necessary to understand that with a 
Maxwellian type of velocity distribution, most molecules have speeds 
close to the average molecular speed. A very small fraction has 
speeds much greater than the average or much less than the average. 

For example, only about two percent of the molecules in a given 
volume of gas have speeds that exceed twice the average molecular 
speed. Thus, it should be apparent that, at high resultant speeds 
(the resultant speed being the component normal to the body U sin 0) , 
the relatively small spread of molecular speeds predicted by the 
Maxwellian distribution become unimportant in comparison with the 
resultant speed. 

This diminishing importance of the spread of molecular velocities 
Is the reason for the constancy of the second term of equation (15) 
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for values of (u/v m ) sin 0 greater than 2. It is interesting to 
note that the term in equation (15) , which accounts for the trans— 
latory molecular motions, has values ranging from 2kT at zero 
velocity to (5/2)kT at values of (U/Vm) sin 6 of 2 and higher. 

This may he compared with the constant value of (3/2)kT given in 
the simplified analysis of Sanger. The increased value at zero 
velocity in the present analysis is due to the fact that the faster- 
moving molecules are "both more likely to collide with the hody and 
also carry more translational energy than do the slower— speed 
molecules. The increase in the term with increasing values of the 
parameter (Tj/Vm) sin 0 is due, as previously inferred, to the 
contribution of the cross-product term of the thermal and mass 
velocities in equation (12), the unintegrated expression for the 
translational energy. The actual difference in the total kinetic 
energy term computed by means of the simplified theory and the more 
exact theory depends upon the relative magnitudes of the terms mO 2 /2 
and tykT. For most of the cases computed in this report, the term mlT /2 
became large in comparison with the term tykT at flight speeds above 
8000 feet per second; consequently, little difference between the 
simplified and more— exact theory would be expected at flight speeds 
greater than 8000 feet per second. 

Physically, It is clear why both the number of molecules 
striking the rear side of the body, and the tranBlatory energy 
incident upon the rear side of the body, decrease with increasing 
values of (u/Vm) sin 0 . The front of the body acts as a shield and 
prevents collisions of the rear side with all molecules except those 
that have absolute velocity components in the direction of flight 
with magnitude equal to the flight speed or higher. Consequently, 
as the flight speed is increased, we should expect that the incident 
energy on the rear of the plate would eventually approach zero. 


Skin Temperatures in the Absence of Solar Radiation 

Figures 7 and 8 show the magnitude of plate tenperatures calcu- 
lated for several altitudes in the absence of solar radiation. The 
marked effect of emitted radiation from the body in lowering the 
plate temperature is striking. For a large range of flight speeds 
and small plate angles, the plate temperatures are actually lower 
than the ambient-air temperature. The curves of temperature of the 
rear side of the body show an eventual decrease with flight speed in 
accordance with the discussion of the preceding paragraph. The effect 
upon body temperature of perfect thermal contact between the front 
and rear sides of a flat plate is shown in figure 9- It will be 
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noted that the plate temperatures are considerably lover than those 
calculated for the front side of an insulated plate, as shown in 
figures 7 and 8. It voulcjl appear that the cooling problem could be 
alleviated by insuring thermal contact between portions of the body 
inclined at positive and negative angles with respect to the flight 
path. Practically, this could be accomplished by boattailing the 
body and using a heavy skin of high thermal conductivity between the 
nose section and the boattailed section. An alternative scheme would 
consist of employing a circulating— type cooling system, abstracting 
heat from the nose section and rejecting this heat at the boattailed 
section. It is also apparent that, from a skin temperature stand- 
point, large positive body angles should be avoided; the body should 
preferably be finely tapered. Fortunately, this dictate coincides 
with aerodynamic requirements. 

The effect of altitude on the temperature of the front side of 
an insulated plate is also indicated in figures 7> 8, and 9- It is 
apparent that the problem of high skin temperatures becomes of negli- 
gible importance above altitudes of 125 miles for flight speeds up 
to 36,000 feet per second and with solar radiation absent. The 
condition of zero solar radiation corresponds, of course, to nocturnal 
flight . 


Internal Cooling 

The effect upon body temperature of internal cooling is shown 
in figure 10 for several flight speeds and for a body ang le of 10°. 
It can be seen that the body temperature decreases rapidly above a 
certain critical value of rate of internal cooling. This effect is 
due to the combined actions of emitted radiant energy and the rate 
of internal cooling. At the lower rates of internal heat removal, 
the emitted radiant energy comprises the major portion of the total 
energy leaving the surface, hence changes in the rate of internal 
cooling have little effect upon the surface temperature. As the 
rate of internal cooling increases, however, and the surface tem- 
perature decreases, internal cooling becomes the controlling factor, 
since the emitted radiant energy is proportional to the fourth 
power of the surface temperature. As the rate of internal cooling 
is increased, a point is finally reached where its magnitude 
approaches the value of the incident molecular energy and when this 
occurs, the surface temperature rapidly approaches zero. 
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The Effect of Solar Radiation, on Skin Temperature 

The variation of plate surface temperatures with flight velocity 
in the presence of solar radiation is shown in figure 11. It can he 
seen that the effect of solar radiation on the surface temperatures 
becomes increasingly important as the flight altitude is increased 
and that the effect is strongest at the lower flight speeds. At 
altitudes greater than 100 miles, the presence of solar radiation 
markedly increases the plate surface temperature and this effect 
becomes practically independent of flight velocity. The implication 
is, of course, that heat transfer by molecular energy exchange Is 
negligible. The strong effect of solar radiation on skin temperature 
at the higher altitudes suggests the desirability of nocturnal flight 
in order to prevent excessive skin temperatures. 


The Effect of Emissivity on Skin Temperature 

The effect of the emissivity of the plate material on the 
surface temperatures is shown in figures 12 and 13 for the condition 
of no solar radiation and in figures 14 and 15 for the condition of 
ma.xi mum solar radiation. It is apparent that it is desirable to 
utilize as high a value of skin emissivity as possible. This conclu- 
sion Is reached because any practical missile must fly over a wide 
range of altitude and speed and, under all conditions of flight speed 
and altitude, a high emissivity tends to reduce skin temperature. 

The effect of emissivity is most marked at the lower altitudes; at 
the higher altitudes the surface temperatures become almost independent 
of the emissivity. This somewhat surprising trend may be explained 
as follows: At altitudes of 100 miles and higher, the magnitude of 
the molecular energy interchange is small In comparison with the 
radiant energy exchange. Thus, neglecting all energy exchange except 
that due to radiation, an energy balance on the plate may be written 
as eJo = eaT 4 . The emissivity e is seen to cancel out and the 
plate temperature is thus dependent only upon the magnitude of the 
solar radiation J c . It should be re-emphasized that these compu- 
tations are valid only for gray-body radiation. 


The Effect of Accommodation Coefficient on Skin Temperature 

The effect of the magnitude of the accommodation coefficient on 
the plate surface temperatures is shown in figures 16 and IT. It may 
be noted that, as the accommodation coefficient is lowered, implying 
increasing difficulty of molecular energy exchange, the plate temper- 
ature likewise decreases. This phenomenon may offer a method of 
cooling, as it may be possible to decrease the accommodation coefficient 
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by altering the surface characteristics of the skin. 


CONCLUSIONS 

The following conclusions can be drawn from the trends shown in 
the preceding analysis: 

1. With solar radiation absent, that is, for the condition of 
nocturnal flight, the skin cooling problem associated with high-speed 
missile flight becomes negligible at altitudes of- 125 miles and higher 
up to steady flight speeds of 36,000 feet per second. 

2. The effect of solar radiation on skin temperatures becomes 
increasingly important as the flight altitude is increased. At an 
altitude of 150 miles and higher, solar radiation is the predomina- 
ting factor that determines skin temperature. 

3. With solar radiation present, the effect of emissivity on 
skin temperatures becomes of decreasing importance as the altitude 

is increased. At altitudes of 125 miles and higher, the skin temper- 
ature at a given flight speed, is practically independent' of the 
emissivity and depends only on- the intensity of the solar radiation. 

4. The strong effect of solar radiation on skin temperatures 
at altitudes above 75 miles suggests the desirability of nocturnal 
flight in order to minimize skin temperatures. 

5. In order to maintain low skin temperatures, the angle of 
Inclination of the body with respect to the flight path should be 
kept as small as possible. This may be done by designing the body 
to be finely tapered and by flying at small angles of attack. 

6. Skin temperatures may be reduced by insuring thermal 
contact between portions of the body inclined at positive and 
negative angles with respect to the flight path. As much surface 
as possible should be inclined at negative angles. Practically, 
this may be accomplished by boattailing the body. 

7. In order for a skin-cooling system to be effective, the 
rate of internal cooling must be of the same order of magnitude or 
greater than the rate at which heat is lost naturally by radiation 
from the body. Internal cooling rates below the natural radiation 
loss rate have little effect on skin temperatures. 
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8 . For a practical missile designed, to fly over a wide range 
of altitudes and speeds# it is desirable to make the emissivity of 
the skin as high as possible. This conclusion# however# is restricted 
to the case of skin surfaces for which the emissivity is independent 
of the wavelength of the emitted and absorbed radiant energy. 

9. A decrease in accommodation coefficient results in an 
accompanying decrease in skin temperature. This phenomenon may offer 
a method of skin cooling if it is possible to lower the accommodation 
coefficient by altering the surface characteristics of the skin. 


Ames Aeronautical Laboratory# 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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Figure 2- Coordinate system used in analysis. 
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Figure 7 The effect of flight velocity and angle of incidence on the 
temperature of the front and rear sides of an uncooled plate. 
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Figure 8 - The effect of flight velocity and angle of incidence on the 
temperature of the front and rear sides of an uncooled plate. 
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Figure 9 - 77ra effect of thermal contact between front and rear 
sides on plate temperature. 
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F igure 10- The variation of plate temperature due to heat removal from the surface. 
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Figure II The effect of solar radiation on the surface temp- 
erature of an uncooled piate. 
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Figure 12 Effect of emissivity on plate temperature assuming 
solar radiation absent. 
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rigure !3 .-Effect of emissivity on plate temperature assuming solar 
radiation absent. 
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Figure 14 - Effect of emissivity on plate temperature assuming solar 
radiation present. 
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Figure 15 rEffect of emissivity on plate temperature assuming 
solar radiation present. 
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Figure !6 -Variation of plate temperature with accommodation coefficient. 
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Figure 17 - Variation of plate temperature with accommodation coefficient. 





